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To show a mechanism leading to the breakdown of a particle picture for the multi-component 
Bose-Einstein condensates(BECs) with a harmonic trap in high dimensions, we investigate the 
corresponding 2-d nonlinear Schrodinger equation (Gross-Pitaevskii equation) with use of a mod- 
ified variational principle. A molecule of two identical Gaussian wavepackets has two degrees of 
freedom(DFs), the separation of center-of-masses and the wavepacket width. Without the inter- 
component interaction(ICI) these DFs show independent regular oscillations with the degenerate 
eigen-frequencies. The inclusion of ICI strongly mixes these DFs, generating a fat mode that breaks 
a particle picture, which however can be recovered by introducing a time-periodic ICI with zero 
average. In case of the molecule of three wavepackets for a three-component BEC, the increase of 
amplitude of ICI yields a transition from regular to chaotic oscillations in the wavepacket breathing. 
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Recently a great number of theoretical and experimen- 
tal efforts have been devoted to Bose-Einstein conden- 
sates (BECs) [J 0. As well as single-component BECs, 
the trapping techniques can create multi-component con- 
densates which involve inter-component nonlinear inter- 
actions. The multi-component BEC, far from being a 
trivial extension of the single-component one, presents 
novel and fundamentally different scenarios for its ground 
state and excitations. In particular, it has been observed 
that BEC can reach an equilibrium state characterized 
by separation of the components in different domains^- 

BEC has a dual aspect of waves and particles. The 
wave nature is high-lighted in the phenomenon of inter- 
ference leading to fringe patterns Q, 0. On the other 
hand, the particle nature of BECs can be seen in typical 
localized states like vortices and solitons. In fact solitons 
were observed in the quasi-one dimensional BEC 0,0. 

Among the works that emphasize a role of the par- 
ticle picture for BEC with a harmonic trap in high- 
dimensions, those of Perez-Garcfa et al. are the most 
noteworthy 0, 0, H, @- We focus on their two important 
assertions. The first one made for a single-component 
BEC is as follows 0: If the phase of BEC wave func- 
tion will be suitably corrected, the center-of-mass for a 
wavepacket displaced from the origin obeys Newtonian 
dynamics and the Ehrenfest theorem is valid even for the 
nonlinear Schrodinger equation (NSE). Furthermore the 
center-of-mass is decoupled from dynamics of the shape 
of a wavepacket. The second assertion is concerned with 
the multi-component BEC0: If the distance between 
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wavepackets with each linked to the individual component 
is much larger than their typical widths, the picture of 
soliton molecules composed of solitonic atoms (wavepack- 
ets) holds for the multi-component BEC with a harmonic 
trap in high-dimensions. The large distance as above is 
guaranteed by a sufficiently large centrifugal force due to 
non-vanishing angular momentum. . 

Perez-Garcfa et a/.'s second assertion is quite interest- 
ing because the solitonic structures like wavepackets are 
believed to be dynamically unstable in two and higher 
spatial dimensions. However, we have several criticisms 
against their assertion: (1) Inter-component interaction 
has a tendency to swell out individual wavepackets and 
breaks their picture of interacting particles, as evidenced 
in Fig^a); (2) Collective coordinates for the width and 
phase of wave packets, which should be coupled with 
the center-of-mass, are not taken into consideration, al- 
though being studied intensively in other works of their 
own[rllFl|. 

In this Letter, choosing the 2-d multi-component BEC 
with a harmonic trap, we develop a refined variational 
principle to derive an effective dynamics for interacting 
wave packets and pinpoint the weak points of the idea 
of Perez-Garcfa et al.^. Then we shall propose a new 
model to restore the picture of soliton molecules in two 
dimensions. In particular, by controlling a parameter we 
shall see a transition of the wavepacket breathing from 
regular to chaotic oscillations. We numerically analyze 
the subject on the basis of a three-component BEC cor- 
responding to the "three-body problem" . 

The multi-component BEC at zero temperature 
is described by NSE (or Gross-Pitaevskii(GP) equa- 
tion). We shall consider a system of n complex fields 
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Fat mode 





Skinny mode 

FIG. 1: (a) Time evolution (up to t = SOlco' 1 ]) of initial 
three Gaussian wavepackets in 3-component BEC with angu- 
lar momentum M = 9 and g — Ag = 2n (computed by using 
Eq.Q). See an outbreak of wave interference; (b) Fat and 
skinny normal modes (schematic). 



tpi(t, r),tp 2 (t,r), . . . , ip n (t, r) ruled by the equations 



are respectively the first- and second-order coefficients of 
Taylor-expansion of the phase w.r.t. r — R(a trivial 
constant phase is suppressed). R, w, a and (3 constitute 
collective coordinates below. 

The expansion of 9 in Eq.J2J is the most natural, al- 
though the existing works E3 employ an expan- 
sion w.r.t. r rather than r — R. The advantage of our 
expansion is that a and (3 turn out to have a transparent 
correspondence with the velocity of center-of-mass and 
the WP breathing, respectively (see JSJ) and © below). 

According to the variational principle, the action func- 
tion deriving Eq. Q is obtained from Lagrangian density 
for held variables: 
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^(M)+X>fe|^(*,r)| 2 V>j(*, 

n (in units of atomic mass m = 1, con- 

— 1 and oscillation period ui^ 1 = 1). 

r" fZ = (x" + y 2 )/2 stands for a harmonic trap and 
Sfe 9jk\4'k(t, r)| 2 is the nonlinear term. The coefficient 
(proportional to the scattering length scaled by con- 
fining length) is tunable by Feshbach resonance. We first 
consider a static case of repulsive nonlinearity, gij > 
and choose gu — g and </y = Ag for i ^ j. 

In the absence of the inter-component interaction(ICI), 
i.e., for Ag = 0, each component has a solitonic state, 
namely Gaussian wavepacket or vortex solution with a 
movable center-of-mass. On the other hand, the non- 
vanishing ICI(Ag 0) among solitons may yield soliton- 
soliton interactions. To reveal a crucial role of the 
breathing nature of above localized structures, we ap- 
ply a refined variational principle and derive from 
the evolution equation for the collective coordinates of 
wavepackets. The trial Gaussian wavepackets(WPs) are 
constructed from the circularly-symmetric ground state 
solution of with ICI suppressed: 
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In fact, the multi-component GP equation is obtained 



from Lagrange equation, 
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We now insert @ into (JSJ) and obtain Lagrangian L 
for the collective coordinates by integrating C over space 
coordinates: 
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Lagrange equations of motion for phase variables a and 
(3 lead to 



exp(«6j(r)), 



= R i 



(8) 



%(r) = a i (r-R i ) + -(r-R j ) T /3 i (r-R i ). (2) 



and 



(One may also choose excited states responsible for vor- 
tices.) R = (X, Y) is the center-of-mass and w is width 
of the circular WP. 



a = (a ,a ) 



(3) 



(9) 



Noting that a. and (3 in (jHJ and 10 (and also all 
higher-order phase coefficients) adiabatically follow R 
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and w, we can rewrite Lagrange equations for the lat- 
ter by eliminating the former: Equation of motion for R, 
A ( JL^.^ 



Jj£- = 0, becomes 



R-t -f- R^ 
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U({B. j },{w j }) = 0. 



(10) 



Similarly for w we have ^ 
reduced to 



(m_) 
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0, which is 
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C7({R i },{ % }) = 0.(11) 



A couple of Eqs. ijlfjl) and (fTT|> indicate: Contrary to 
Perez-Garcfa et a/.'s theory that ignored the role of WP 
width0], the width and center-of-mass are strongly corre- 
lated, leading to a breakdown of the picture of a solitonic 
molecule based on the center-of-masses alone. To exam- 
ine the problem in detail, we consider the two-component 
BEC in two dimensions and explore the fate of a WP 
breathing against ICI (Ag). Effective Lagrangian leading 
to Eas. ((TUjl and (|l"l"1) is now expressed as 
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Let us define the center-of-mass of two components 

R.2-R.1 



R = Rl +fi 2 and the relative displacement Q 

V 2 



V2 



and suppress the global translational degree of freedom 
(Ro = Ro = 0). Owing to rotational symmetry the 
angular momentum M is a constant of motion: M = 
dL "f f — Q 2 9 with 9 as a polar angle for Q. Choosing 

ou 

a synchronous width dynamics (w\ = W2 — w) with the 
canonical momentum p w — -^(L e ff — M6), the effective 
Hamiltonian is given by 

Has = \PQ 2 + \pw+V{Q,w) 



with 



V(Q,w) 
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Equation 1|13|) indicates: In the absence of ICI(A = 0), 
w and Q show independent regular oscillations around 
the minima of individual potentials, w = wq = (1 + 
g/2n) 1 / 4 and Q = Qq = y/M, respectively. Interestingly, 
however, for any set of values M and g the typical os- 
cillation frequencies ujq for Q and oj^ for w are always 
degenerate: uiq — w° = 2 (in unit of frequency of the 
harmonic potential, u>). Therefore, by switching on a 
small but non-vanishing ICI, dynamics of w and that of 



Q will be strongly mixed to remove the degeneracy and 
normal modes, fat (optical) and skinny (acoustic), are 
formed. In the fat mode, Q and w change in anti-phase, 
namely the decrease of the relative displacement is ac- 
companied by the blowup of WPs, while in the skinny 
mode, Q and w change in phase (see Fig^b)). The 
emergence of the fat mode brings about the wave in- 
terference between adjacent WPs as shown in Fig^a). 
This provides a mechanism leading to the breakdown of 
Perez-Garcfa et aZ.'s picture of interacting particles. 
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FIG. 2: Upper panel: Time evolution of WP width (wj). 
g = 2n, Ag — 3 and angular momentum M = 9. Q is a tunable 
parameter, (a) f2 = 0(dotted line) and Q — 100(solid line), 
(b)f2 = 3,(c)f2 = 10; Lower panel: Time evolution (up to 
t = 100[w -1 ]) of initial three Gaussian WPs in 3-component 
BEC with g = 2tt, Ag = 3, Q = 100 and M = 9. See the 
persistence of a particle picture. 

To suppress the wave interference between adjacent 
WPs and recover the particle picture, we proceed to in- 
troduce a time-periodic ICI with zero time-average and 
the amplitude Ag, 



gij = Agcos(fit) [i^j) 



(14) 



rather than sticking to the static ICI. Here the frequency 
scale w° (= 2) is essential. If O 3> the width will feel 
actually no ICI, showing a stable and regular breathing 
proper to a single component BEC with a harmonic trap. 
For n near but larger than , one can expect a suffi- 
ciently long stable breathing until the time 27r/w°) 
that a blowup of WP will start. For fl equal or less than 
uj^, the breathing oscillation will quickly blow up leading 
to the wave interference. 

To verify the above conjecture, we have numerically 
integrated 3-component NSE(GP equation) with a har- 
monic trap in two dimensions in Q with use of split- 
step Crank-Nicholson method 01 • Initial data are co- 
ordinates Xj = Rcos(2ir(j — 1)/3),Y} = i?sin(27r(j — 
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FIG. 3: Time evolution of WP width (wj) (upper panel) and 
corresponding power spectrum (lower panel), g = 2ir, SI = 5 
and M = 9. A$ is a tunable parameter. (a)(a')A<; = 
0,(b)(b')A<? = 2, (c)(c')A 5 = 4. 



l)/3) and velocities Vjx = —Ruism(2ir(j — l)/3), VjY = 
Raj cos(27r(j - l)/3) with R 2 ui = M/3 (M: angular mo- 
mentum) and j = 1,2,3. In short these are placed on 
the vertices of an equilateral triangle. We have also cal- 
culated expectation values of the widths using w? x = 
J((x- < x >f\^\ 2 )dv,w 2 Y = J((y- < y >) 2 \^\ 2 )dr 
for different times. In Figs. |21 and [3] the mean value 
Wj = (wjx + wjy)/2 is plotted as a function of time. 
Note: even for a single-component NSE^, Q, the ini- 
tial circular symmetry of WP is slightly broken during 
the time evolution, due to the nonlinear intra- component 
interaction; Each of the expectations Wjx and Wjy is de- 
scribed by a superposition of modes with nearly-equal fre- 
quencies w° and + g/4n)/(l + g/2n)uj^, and shows 
a small beating, which however can be suppressed by 
choosing the above mean value Wj . 

Figure[2Ia) shows a blowup of the breathing oscillation 
for the static ICI (SI = 0) and a recovery of stable oscilla- 
tion for SI ^> oj^, and Figs^b) and|2^c) show oscillations 
with no blowup for ft > w°. For a fixed value of SI in 
the range SI ~ w", we see a transition from regular to 
chaotic breathing oscillations when the amplitude Ag is 
increased. In fact, in the case of SI = 5, the oscillation is 
irregular for Kg = 2 and 4 (see FigsOJb) and|3fc)), while 
it is regular for Ag — (see FigEfa)). The correspond- 
ing power spectra show a transition from a distinguished 
line structure in FigEHV) to broad ones in FigEJb') and 
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FIG. 4: Phase diagram in Ag-Q space, (g = 2ir and M = 9) 
x,0,n stands for blowup, chaotic and regular regions, re- 
spectively. (Data are increased in the left half of the diagram 
for use of detailed discussions.) 



He'). 

We have systematically investigated the wave dynam- 
ics starting commonly from the equilateral triangle con- 
figuration as above. Figure0]is a phase diagram in Ag-fl 
space, which shows three distinct regions, (i) blowup, (ii) 
regular without blowup and (iii) chaotic without blowup. 
We find the decrease of the blowup region as SI is in- 
creased beyond w° and the emergence of chaos without 
the blowup as Ag is increased. Contrary to our simple 
expectations, the phase diagram involves rich informa- 
tion to be explored: The regular and chaotic oscillations 
without blowup are observed even for SI < oj^, and the 
blowup occurs at resonances, SI = and 2w° , even for 
a small amplitude Ag. However, Fig0] clearly conveys 
that the three-component BEC with the oscillating ICI 
restores a picture of interacting particles and can exhibit 
the transition from regular to chaotic motions without 
the blowup. Thus we have obtained a model of conser- 
vative chaos constructed from the three-component BEC 
with a harmonic trap in 2 dimensions. 

To develop Perez-Garcia et al.'s idea, we have exam- 
ined the multi-component repulsive BEC in a 2-d har- 
monic trap. In the absence of the inter-component in- 
teraction (ICI), the wavepacket(WP) breathing and the 
motion of relative distance between WPs have oscilla- 
tion frequencies always degenerate for any set of the an- 
gular momentum and intra-component interaction. The 
non-zero ICI removes the degeneracy and induces the fat 
mode that breaks a picture of soliton molecules. We have 
therefore proposed a new model with a time-periodic ICI 
with zero time-average, which can elongate the time in- 
terval without any blowup of WPs. Provided that the 
frequency is near the characteristic breathing frequency 
(oj^,), a transition from regular to chaotic oscillations oc- 
curs as the amplitude of ICI is increased. For attrac- 
tive BECs without a trap or in free space, the oscillating 
nonlinearity with non-zero average is known to stabilize 
high-dimensional WPs[l3L FLU]. In a similar way, in the 
case of the multi-component repulsive BEC with a trap, 
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the analogous picture holds for high-dimensional WPs 
mutually interacting through time-periodic ICI with zero 
average. The amplitude of ICI is found to control the 
chaoticity of breathing motions of WPs. The analysis 
of chaos in space coordinates and the comparison of the 



present results based on NSE with the effective dynamics 
with a few degrees-of-frecdom in Eqs. IjlUI) and (|llf> will 
constitute subjects which we intend to study in future. 

K.N. is grateful to M. Lakshmanan for valuable con- 
versations in the early stage of the present work. 
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